Dynamic process of free space excitation of asymmetry resonant microcavity 
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The underlying physics and detailed dynamical processes of the free space beam excitation to 
the asymmetry resonant microcavity are studied numerically. Taking the well-studied quadrupole 
deformed microcavity as an example, we use a Gaussian beam to excite the high-Q mode. The 
simulation provides a powerful platform to study the underlying physics. The transmission spectrum 
and intracavity energy can be obtained directly. Irregular transmission spectrum was observed, 
showing asymmetric Fano-type lineshapes which could be attributed to interference between the 
different light paths. Then excitation efficiencies about the aim distance of the incident Gaussian 
beam and the rotation angle of the cavity were studied, showing great consistence with the reversal 
of emission efficiencies. By projecting the position dependent excitation efficiency to the phase 
space, the correspondence between the excitation and emission was demonstrated. In addition, we 
compared the Husimi distributions of the excitation processes and provided more direct evidences 
of the dynamical tunneling process in the excitation process. 
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I. INTRODUCTION 

Optical microcavity has drawn many attentions owe 
to its high intracavity field intensity [1]. Light-matter 
interactions, such as laser (H, nonlinear optics [3], cav- 
ity quantum electrodynamics 0, (H, sensors |6i] and op- 
tomechanics Q, can benefit from the microcavities. No 
matter these applications are based on active or passive 
microcavities, energy exchange with the outside is essen- 
tial. Variously shaped asymmetry resonant microcavities 
(ARCs) (quadrupole spiral stadiu m I 111 , Li- 

macon [l2f, and half-quadrupole-half-circle [l3|,|l3|) were 
extensively studied recently since light can be collected 
more efficiently and conveniently through free space in 
those cavities than in symmetry circular or spherical cav- 
ity. 

However, most studies put emphasis on the directional 
emission for efficient collection (HI, lol-ITBI , but the inverse 
processes, i.e. the free space excitation [T], have not 
been fully explored. An's group performed quadruple 
droplet pumping experiments, then put mode-coupling 
theory and numerical simulation forward to explain the 
experimental phenomenon p!6l-[T8| . Their studies focused 
on the laser emission in a full-chaotic cavity by a plane 
wave excitation, and the process of chaotic transport in 
the pumping process was studied indirectly from the laser 
threshold. However, threshold also depends on the over- 
lap between the pumping mode and lasering mode, and 
there is also mode competition in the multimode cavity. 
The laser threshold may not linearly depend on the exci- 
tation efficiency of pumping mode. Therefore, the direct 
evidence of the efficient pump to a passive high-Q modes, 
and the studies of underlying principle are required. In 
addition, more perspectives of the free space coupling. 
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such as the focused beam pump and the spectrum, are 
lacking. 

In this paper, the free space excitation of an ARC by a 
focused Gaussian beam is studied numerically by bound- 
ary elements method (BEM). Then the field distribution, 
mean intracavity intensity (Mil) and the transmission 
spectrum varying with frequency of the incident beam 
are also studied. Fano-type lineshape of spectrum is in- 
vestigated from interference point of view. Excitation 
efficiency represented by Mil is studied by adjusting rel- 
ative position of the incident beam and the cavity. A 
comparison between excitation and emission is performed 
with the help of ray optics. Finally, resorting the Husimi 
projection of the wave field and the projection of the Mil 
distribution, we discuss the physics mechanism of the free 
space excitation in several cavities with different degree 
of deformations. 



II. THE BASIC DYNAMICS OF QUADRUPOLE 
SHAPED MICROCAVITY 

Here, we focus on the two dimensional (2D) quadrupole 
microcavity, which has been studied extensively in the- 
ory and experiments [To']. Basic features and underly- 
ing physics found in this specific cavity can also extend 
to other ARCs, like the half-quadrupole-half-circle mi- 
crocavity. The boundary of quadrupole microcavity in 
polar coordinate is described by r(^) = R{1 + £ cos 20) 
with radius R and degree of deformation e In the 
deformed cavity, the ray dynamics shows quasi-periodic, 
periodic and chaotic trajectories, corresponding to the 
regular KAM torus, islands, and chaotic sea in the phase 
space presentation. Fig. [TJc) (black dots) shows the ray 
dynamics in quadrupole cavity with deformation £ = 0.1 
as an example. Energy of a high quality factor mode is al- 
ways stored in the vicinity of the regular structure above 
the critical line by total internal reflection [Fig. HJa)], 
part of energy tunnels to the nearby chaotic zone and 
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Figure 1. (a) Period-4 rectangle obits. The ray trajecto- 
ries are calculated from the location of maximal points in 
the Husimi projection in (c). Considering the Goos-Hanchen 
shift, clockwise (red solid line) and counter clockwise (blue 
dash line) obits split with each other, and the obits are not 
closed, (b) Near field intensity of the order-2 mode, (c) 
Poincare Surface of Section (SOS) of the quadrupole cavity 
(black dots) and Husimi projection (green areas) of the WGM 
illustrated in (b). In the SOS, ray is recorded as a point with 
coordinate (0, sinx), where is the polar angle of incident 
point, and x is incident angle on the boundary, respectively. 
The whole SOS of a quadrupole cavity is mixed with Kol- 
mogorov-Arnold-Moser (KAM) torus (curves near the up- 
per and lower boundary of SOS), islands (close curves), and 
chaotic sea (collection of other dispersed points) structure. 



finally leaking out the cavity by reflections. The reflec- 
tions condense in a few trails, which are called manifold, 
in phase space [Fig. 5(a)] [20]. The asymmetrical nature 
of the manifold results in the directional emission. 

Ray model just gives basic physics insights of the 
ARCs. For more underlying physics and precise details, 
we should also consider the aspect of wave nature of light. 
The wave functions of high-Q WGMs are calculated nu- 
merically (solving Maxwell function using boundary ele- 
ments method (BEM) [2T|). Setting uniform cavity re- 
fractive index n = 1.45, we search the resonances with 
quick root searching method [22] at kR ~ 50, where kR 
is also the dimensionless cavity size scale to beam wave- 
length. Figure [T](b) is a eigenmode of quadrupole cavity 
near kR = 50, it is known as the order-2 mode for exist- 
ing two intensity peaks along the radial direction. We can 
see the ray-wave correspondence by Husimi projection of 
the wave function in phase space, as is shown by the 
green contour in Fig. 1(c). Eight obviously condensed 
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Figure 2. A schematic diagram of an incidence beam, a cav- 
ity, and a detector used in the simulation. Note that this 
configuration is conjugated with former experiments as aim 
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regions are clockwise and counter-clockwise period-4 is- 
lands, which corresponds to the period-4 orbits in Fig. 
1(a). Noting that the orbits shows a displacement at 
each reflection, which is attributed to the Goos-Hanchen 
shift [23|. In the following studies, we focus on this mode 
to study free space excitation. 



III. THE CONFIGURATION OF GAUSSIAN 
BEAM EXCITATION 

Now, we turn to study the free space coupling to the 
WGM [Fig. 1(b)] by Gaussian beams [Fig. [g. From the 
perspective of wave optics, incident focused 2D Gaussian 
beam reads UM 



A(x,z) 



1 



exp[- 



■ ikz] 
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Here ( = 2z/b is the dimensionless longitudinal location 
defined in terms of the confocal length b = kuo^^ and Co'o is 
the beam waist. This Gaussian beam is set as cjq = 0.3i?, 
and focus on 2; = line. The relatively narrow waist of 
the Gaussian beam can be obtained in an experiment by 
a convergent lens, and benefit the coupling efficiency. 

The coupling efficiency depends on the excitation con- 
dition of the Gaussian beam, i.e. the direction of the 
beam and the excitation point on the cavity boundary. 
So, we rotate the cavity around the original point, and 
shift the incident beam parallel along x axis. Thus, the 
coupling config can be represented by two parameters 
iO^x)^ where is angle from axis z to major axis of the 
quadrupole, and x is the aiming distant between the ax- 
ial line of the Gaussian beam and the central point of the 
cavity, which is just the numerical value of x coordinate 
of the axial line in our coordinates system. 

In the first place, an on-resonance and an off-resonance 
excited patterns of light intensity in cavity with £ = 0.1 
are shown in Figs. [3ja) and[3l^b). All wave simulations 
are performed by BEM, which solves the Helmholtz equa- 
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Figure 3. False color intensity patterns of excited cavity using 
on-resonance (a) and off-resonance (b) incident beam. White 
lines are the boundaries of the cavity, (c) Angular distribu- 
tions of outward energy flow through a circle with radius 1.8R 
correspond to the pattern (a) (red line) and (b) (blue line), 
respectively. A positive value indicates the outflow, a nega- 
tive value indicates inflow. Parameters: kR = 49.726 (on), 
49.9 (off), = 45°, X = 0.97, e = 0.1. 



tion of the electromagnetic field with particular bound- 
ary conditions. For simplicity but without loss of gen- 
erality, only transverse magnetic mode is taken into ac- 
count. In this polarization, electric field vector is vertical 
to the cavity plane. The field amplitude and its deriva- 
tive are both continuous on the interfaces of dielectrics. 
Incident beam frequency kR = 49.726 (49.9) is chosen 
for on-resonant (off-resonant) coupling, = 45°, and 
X = 0.97 are set for getting high intracavity intensity. 
As is shown in Fig. [3l normalized 2D Gaussian incident 
beam is partly coupled into the cavity. 

In the on-resonance case, notable intensity lies near 
the cavity boundary forming the whispering gallery like 
mode, which is consistent with a traveling wave mode 
corresponding to the order-2 resonance mode at the same 
kR [Fig. [Hb)]. There is no interference fringe along the 
(/) in Fig. 3(a), which is due to only counter-clockwise 
traveling WGM excited in such excitation configure with 
a focus Gaussian beam pumping away from the center of 
the cavity. When it comes to off-resonance case, the cav- 
ity field is very different from the eigenmode field. The 
incidence beam will be refracted into the cavity, reflect 
along the cavity boundary several times, and finally leak 
out by refraction [25]. Due to the chaotic ray dynamics, 
the rays show large expansion in the r direction. This 
large expansion can also be interpreted as excitation of 
multiple low-Q modes. The field of which shows higher 
radial quantum number. These multiple interference also 
occurs in the on-resonance case, but is not visible due 
to high resonant field intensity. Note that the incident 
Gaussian beam is also not visible in the on-resonance 
condition, which is due to the greatly enhanced cavity 
field making the relatively low intensity of pumping field 
invisible. In order to evaluate the intensity in the cavity. 



we introduce the Mil as average intensity index. When 
the maximum intensity of the incident Gaussian beam on 
the waist line is set to a unit as reference intensity, the 
Mil is proportional to coupling efficiency and the Q of a 
mode. For a particular mode, such as order-2 mode that 
we concerned, Q is fixed, so the Mil represents the cou- 
pling efficiency directly. The Mils are 5.02 and 0.23 for 
on-resonance and off-resonance excitation, respectively. 
Mil of on-resonance excitation is much larger than that 
of off-resonance excitation. So in the figure of relative 
intensity pattern [Fig. [3](a)] Gaussian beam leaves only 
a dim trace against to the strong intracavity intensity. 

Moreover, angular distributions of radial energy flows 
of on-resonance and off-resonance [Fig. ^c)] were calcu- 
lated by the Poyntin^ vector in radial direction p ■ ft = 
Im(£^ * ■dvE)/2ujjii j2a| on a circle with radius of l.SR 
and center on the original point, where dy = iT • V is 
the normal derivative on radial direction n, uo is angu- 
lar frequency, and /i is relative permeability. The flow 
is negative at 30 degree, which shows the position of 
the incident Gaussian beam, while the other region is 
always positive which indicates the outgoing field. For 
both on- and off- resonance, the net flow as the sum of 
energy flow in all direction is zero, which indicates that 
the energy conservation is preserved, i.e. the ideal pas- 
sive cavity (with real refractive index) does not produce 
or consume energy. It can be seen from the Fig. ^h) 
that the energy flow distribution of off-resonance excita- 
tion mainly embodies the redistribution of the incident 
energy by cavity-caused scattering. In the on-resonance 
case, there is also the emission from resonant modes. So, 
two emission peaks of the order-2 CCW mode emerge 
near 60° and 240° on detector circle in Fig. [D^c), respec- 
tively. The interference between the emission of high-Q 
mode and reflection light of multiple beam also give rise 
to the modulation of distribution of energy flow in other 
directions. 



IV. THE TRANSMISSION SPECTRUM 

Similar with experiment, through a frequency scan to 
the incident beam, we can obtain the spectrum by using 
a detector with fixed relative positions of incident Gaus- 
sian beam and microcavity [with = 45° and x = 0.97]. 
The energy [Fig. [J] red curves] is the integration of the 
outward energy flow through the line detector. In addi- 
tion, the Mil spectrum was also recorded [Fig. S] blue 
curve]. The transmission spectrum near kR = 50 shows 
the shape of the periodic low-frequency harmonized wave 
superimposed with high-frequency spikes. Therein, the 
low-frequency harmonized modulation is interference of 
light that undergo different quasi-period rounds ^5|. The 
modulation period is about 0.7 1/n, i.e. roughly corre- 
sponding to an optical distance of one round trip. 

On the periodically modulated spectra background, 
there are many sharp symmetry dips, peaks and asym- 
metric Fano-type line shapes in the transmission spec- 
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Figure 4. (color online) Spectrum of energy flow (red curves) and mean intracavity intensity (blue curves). The enlarged scales 
are four typical lineshapes in transmission spectrum. 



trum. Comparing with the Mil spectra, these sharp 
variations correspond to a magnified spikes in MIL The 
correspondence indicates that the formation of those line 
shapes is related to the on-resonance excitation. We note 
that the Mil spectrum clearly shows two sets of resonance 
peaks. The set with larger intervals of adjacent lines 
and larger Mil is order- 1 modes and the other is order-2 
modes which are concerned in this paper. Different from 
the Lorentz line shape in the transmission spectrum of 
waveguide coupled cavity, here the line shapes are formed 
by the interference between emission light of resonance 
CCW mode and the scattered light by multiple beam 
interference. 

There are two different paths for light to transmit from 
source to detector: (1) the light coupling to the high- 
Q cavity mode and emitted from the mode, i.e. direct 
scattering from high-Q mode (2) the light transmitted 
after multiple beam interference, or we can say scattering 
from multiple low-Q modes. The final collected energy is 
given by 



E{5) = Emuie'' 



E e'^ ^ 

-'-^res^ . c- / 

10 1 - 



1 



(2) 



Here, 5 is the frequency detuning, n is the linewidth of 
high-Q mode, a and {3 are phase factors that denotes 
the phase variance between source and detector, Emui 
and Eres are the amplitude (real number) of the multi- 
ple beam interference and resonance scattered light. For 
the multiple beam interference, the low-Q modes involved 
have linewidth much larger than a^, thus the phase a is al- 
most constant ior 5 / hi 1. From the formula, the phase 
difference A = a — /3 determines collected transmission 
spectra, (i) A = gives symmetry Lorentz lineshape dip 
(the 4th line of the upper panel in Fig. |4]). (ii) < A < tt 
is left-low-right-high asymmetry Fano shape (the 3rd line 
of the upper panel in Fig. |4]). (iii) A ~ tt gives symmetry 
Lorentz lineshape peak (the 2nd line of the upper panel in 



Fig. Hj). (iv) TT < A < 27r shows left-high-right-low asym- 
metry Fano shape (the 1st line of the upper panel in Fig. 
Hj). The A varies with the frequency of the incident light, 
thus we observed all four different lineshape types in the 
transmission spectra. Then, interference line shapes are 
formed. But unlike the waveguide coupled case, the en- 
ergy flow irradiates into all directions in free space. So 
the spectrum collected by a detector in a particular di- 
rection has no energy conservation relationship with Mil, 
which makes the evaluation of pumping efficiency from 
transmission spectrum invalid. That is why we turn to 
choose the Mil as pumping efficiency indicator. In addi- 
tion, the phase also varies according to different position 
of the detector making different lineshapes (27j . 



V. DETAILS OF EXCITATION PROCESS 

The resonant mode of deformed quadrupole microcav- 
ity has its favorite emission direction. Hence the free 
space coupling should be the most effective under a par- 
ticular pumping condition, as a reversal of directional 
emission from deformed microcavity [27]. Here, we stud- 
ied the pumping with different parameters (6>,x) [Fig. 2], 
and recorded the Mils of four quadrupole microcavities 
(e = 0.00,0.05,0.10,0.15) as illustrated in Fig. Mils 
of circular cavity (e = 0) do not change with the rota- 
tion angles because of the rotational symmetry of the 
cavity. Moreover, the maximum of Mil occurs beyond 
the boundary, which indicates the most effective way of 
excitation the circular cavity is barrier tunneling. The 
slightly deformed quadrupole cavity with deformation of 
0.05 has angle selectivity. Nevertheless, the main excita- 
tion manner is still tunneling, since the maximum pump 
position is outside the cavity too. In addition. Mil dis- 
tribution of this cavity shows a roughly horizontal sym- 
metry, which is a result of the tunneling coupling show- 
ing CW and CCW symmetry. However, in quadrupole 
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Figure 5. Pump efficiency distributions of quadrupole cavi- 
ties with e — 0.00,0.05,0.10,0.15 in real space position co- 
ordinates The frequency of incident beam fixed on 
resonance with the q = 2 mode. White dash fines are the pro- 
jection boundary of cavity in x axis. Because of the decrete 
symmetry, it is sufficient only 0° < ^ < 90° is investigated. 
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Figure 6. (a) Emissions of the quadrupole cavity with s — 
0.1 predicted with ray optics. The points in the figure are 
rays whose refiection angles are larger than the critical total 
refiection angle. So, collections of refractions of those rays 
constitute the emission picture, (b) Phase space presentation 
of s = 0.1 panel of Fig. [5l 



cavity with deformation of 0.10, the highest Mil area is 
located near the boundary projection line, i.e. incident 
beam is tangential to the boundary of the cavity. Be- 
cause the width of the incident beam is not zero, half of 
the lights enter into the cavity by refraction. The size of 
the strong excitation region in the 6 — x figure shrinks 
significantly comparing with the slightly deformed case. 
This localization suggests strengthening of selectivity of 
pumping angle. In addition, the Mil of highest excitation 
region in the upper part of the figure is much stronger 
than that of the second highest excitation region lying 
in the lower part. Different from the barrier tunneling 
coupling where the tunneling rate only depends on the 
local curvature and is same for CW and CCW wave, the 
dynamical tunneling coupling is sensitive to the energy 
distribution of the incident beam in the chaotic sea of 
phase space. The dynamics of rays incident in these two 
regions are non-identical due to the asymmetric unstable 
manifold in the chaotic sea, thus the energy distributions 
in cavity for these two excitation conditions are different, 
and finally raises different excitation efficiency to high-Q 
modes by dynamical tunneling. By increasing the defor- 
mation to 0.15, the highest Mil zone is mainly located 
inside the boundary projection line, which suggests re- 
fraction excitation. This is the only optimal position as 
the dynamical tunneling dominats the pumping process. 
And the Mil is visible when a beam shines on the cavity 
which is due to low contrast to the optimal coupling. The 
Mil is always nonzero when light shines on the cavity. 

As the inverse process of emission, excitation of an 
ARC should show some correspondence with directional 
emission. For example, comparing our pumping effi- 
ciency with the emission figure obtained in Ref . [15| , they 
show great consistence. In order to express this correla- 
tion directly. Mil distribution [e = 0.1 panel of Fig. [5|) 



is projected to phase space [Fig. [6jb)] by ray trajectory. 
First, each point in real space distribution is treated as 
rays. The rays impacted on the cavity (points between 
the white dash line) can refract into the cavity. Then, 
record pairs of refracting positions and the sine of re- 
fraction angles (0, sin x) with pump efficiency denotative 
intensity. For comparison, the phase space presentation 
of emission property of the same cavity is obtained by 
ray simulation [Fig. El^a)]. Those two distributions in 
phase space show good consistence. Thus, the excitation 
is verified numerically as inverse process of emission. 

To investigate the detail of the excitation process and 
get a deep physical insight of the free space coupling, 
Husimi projections of the excitation (purple background 
picture in Fig. [7j) are calculated for different pump con- 
dition. As demonstrated in Fig. [3g), Husimi of the opti- 
mum excitation is consistent with the Husimi of order-2 
mode [see Fig. [TJc) green dots] in CCW half phase plane. 
It shows that the incident energy is coupled into the cav- 
ity and forms resonant traveling mode [Fig. [71(a)]. The 
excitation process is also demonstrated in the case of de- 
tuning light {kR = 50) with the same geometry [Fig. 
[JKb)]. The Husimi projection of the detuning excitation 
has four dense points which have good superposition area 
with the resonance islands. Next, the excitation field cor- 
responding to the second maximum Mil configuration, 
i.e. the optimal excitation on the other side, is obtained 
[Fig. EJc)]. It is a resonant order-2 CW mode as is shown 
in SOS [Fig. [IJc)]. The right side patterns in Fig. [71 are 
the distributions of excitation CCW mode by tangential 
incident beam at the rotating angle of 30° [Figs. [T^d), 
(j)], 20° [Figs. [TKe), (k)] and 10° [Figs. [T^f), (1)], re- 
spectively. The corresponding Husimi projections refiect 
the change of coupling strengthens between the incident 
beam and the resonance mode from strong to weak. The 
mode built in 30° case is merely strong, while the mode of 




Figure 7. False color intensity patterns of excited cavity ((a), (b), (c), (d), (e), and (f)) and corresponding Husimi projection 
((g), (h), (i), (j), (k), and (1)). The aim distances and the rotation angles are labeled with corresponding letters in e = 0.1 
panel of Fig. O The coordinates of real space patterns and phase space projection are the same with Fig[3l^a) and Fig. Hl^c), 
respectively. 



20° only has four weak points in the SOS. As for the 10°, 
only very subtle mode energy emerges in the SOS. As we 
can see in the phase space, the more overlap between the 
incident beam and the resonant mode there is, the higher 
excitation Mil will be established. The overlap owes to 
the wave character of the beam. Despite the penetration 
into islands is forbidden for the refraction rays from the 
outside, the dynamical tunneling effect allows it. 



space. Mil distributions of cavity with different defor- 
mations also exhibit two main excitation processes, i.e. 
barrier tunneling and dynamic tunneling. Besides the 
frequency matching, the sufficient overlap between the 
incident beam and the resonance mode is the key factor 
for efficient excitation. 
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VI. CONCLUSION 

We have gotten the inner and outer field of the exci- 
tation cavity by BEM. Scattering energy flow was also 
calculated with the Poynting vector. The transmission 
spectrum of a certain direction can have asymmetric line- 
shape because of the interference. Mil distribution agrees 
well with the emission picture in reverse manner in phase 
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